P- and T-violating ttNN form factor 



S. Kondratyuk and P. G. Blunden 

Department of Physics and Astronomy, 
University of Manitoba, Winnipeg, MB, Canada R3T 2N2 
(Dated: February 9, 2008) 

Abstract 

The form factor of the parity and time-reversal violating (PTV) pion-nucleon interaction is cal- 
culated from one- loop vertex diagrams. The degrees of freedom included in the effective lagrangian 
are nucleons, pions, r], p and lv mesons. We show that by studying the form factor one can constrain 
the PTV meson-nucleon coupling constants. We evaluate the mean square radius associated with 
the PTV irNN vertex. Using the mean square radius, we estimate the effect of the PTV ttNN 
vertex on the neutron electric dipole moment, and find a very small correction. We also extract 
the renormalisation group (3 function and use it to discuss evolution of the PTV ttNN coupling 
constant beyond the hadronic mass scale. 

PACS numbers: 11.30.Er, 13.75.Gx, 24.80.+y 
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I. INTRODUCTION 



The parity (P) and time-reversal (T) nonconserving pion-nucleon couphng is a theoretical 
ingredient in many calculations of electric dipole moments (EDMs) of various systems, such 
as the neutron, deuteron, various atoms and molecules The continuing interest in the 
problem of EDMs, and its relevance to possible extensions of the Standard Model 0], is 
a major motivation for a deeper understanding of the P- and T-violating (PTV) nNN 
interaction. Usually in calculations one uses a point-like PTV ttNN interaction, neglecting 
a form factor associated with effects of nonlocality. 

In this paper we calculate such a form factor from meson loop corrections to a bare PTV 
TiNN vertex. The loop integrals are evaluated in a relativistic field theoretical approach 
based on an effective lagrangian comprising the nucleon, pion, p, uj and rj mesons. The 
lagrangian consist of a strong interaction piece involving in principle well-known coupling 
constants, and a PTV piece involving unknown coupling constants. Our aim is to study 
various features of the calculated PTV ttNN form factor - for example, its dependence on 
the square of the momentum transferred by the pion to the nucleon, the associated mean 
square radius, etc. Using renormalisation group methods in a one-loop approximation, we 
also consider possible scenarios of the high-energy evolution of the PTV ttNN constant. 
The latter problem in particular is motivated by the puzzle of the extremely small value 
of the charge and parity (CP) violating 6 term in the QCD lagrangian [5]. Throughout 
the paper, we do not assign specific numerical values to the unknown PTV meson-nucleon 
coupling constants. We do, however, make general assumptions as to the relations among 
them (e. g. choosing the constants to be of comparable magnitude or assuming some of them 
to be suppressed) and analyse consequences of these assumptions for the PTV ttNN form 
factor, the mean square radius and the high-energy scale evolution of the PTV coupling 
constant. 

Our results could be used to study effects of the PTV ttNN form factor in calculations 
of electric dipole moments of the nucleon. We will make a tentative estimate indicating that 
the PTV irNN form factor introduces a very small correction to a leading chiral contribution 
to the neutron EDM. To our knowledge, a measurement of a PTV irNN form factor has 
not been attempted so far. We will show that such a measurement could, in particular, put 
some constraints among the PTV meson-nucleon coupling constants. Since at present such 
experimental applications are rather speculative, in this paper we will focus on the PTV 
form factor itself, leaving specific calculations of related measurable PTV amplitudes for 
future research. 



II. EFFECTIVE LAGRANGIAN 

The interaction among hadrons in our model is described by the sum of the strong 
lagrangian Cs and the P-, T-violating lagrangian £fy [§]: 

C = Cs + (1) 

where 
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+ {o^^l.d. ^fr . p, + gfpl + ^(2)r3p2] } ^ 

+ {a^^755. + £^^3^J } iV (3) 

(our constants g^^^ and ^^^j, are equal to — ^j^^j, and 3^^^),, respectively, as defined in 
Ref. jsj). The nucleon field is denoted as A^; the pion, p, oj and 77 meson fields as vr*, p^, 
cj^ and r], respectively, with the lower indices denoting Lorentz vector components and the 
upper indices isovector components; the bold font indicates isovectors and Xj are the Pauli 
isospin matrices.^ The hadron masses are: M = 0.938 GeV for the nucleon, m = 0.138 GeV 
for the pion, nip = 0.77 GeV, = 0.78 GeV, = 0.55 GeV for the other mesons. 
The coupling constants in the strong lagrangian C are chosen as in Ref. g^ = 13.07, 
g^i = 2.24, gp = 2.75, Hp = 3.7, g^^ = 8.25, = —0.12 and ^fp^ = 6.07 (from the width of the 
p resonance ^). 

Retaining all isospin structures in the PTV ttNN (and other meson-nucleon) vertices 
would involve many unknown coupling constants, which would make results of loop calcula- 
tions based on the lagrangian Eq. difficult to interpret. Therefore, to keep our analysis 
tractable, in the following we will consider only the ~ g^^ isospin structure in the PTV 
TcNN vertex, which involves a charged pion and therefore appears in the calculations of the 
neutron EDM j^, [6|. We will use the simplified notation for the renormalised (physical) 
coupling constants: = g^°\ Cp = gf\ = gl^\ = gl^\ Note that the PTV vrA^A^ 
coupling constant in the lagrangian Eq. ([3]) is a bare constant which will serve to cancel 
infinities of the loop integrals in the course of renormalisation described in the next section. 

We use the pseudovector (pv) strong ttNN vertex in Eq. ([2]). The pseudoscalar (ps) 
structure igj^Nj^T ■ ttN is also used in the literature on effective lagrangian calculations of 
effects of parity and time reversal violation [1, |6|. The pv structure of the nNN vertex was 
used in the recent chiral calculation [3] of the nucleon electric dipole form factor. While the 
ps and pv couplings are equivalent on shell, the pseudovector is more consistent with chiral 
constraints, which motivates our choice. In the course of the discussion below we will specify 
important differences and similarities between the results obtained using the pv and the ps 
strong nNN vertices. 



III. ONE-LOOP VERTEX DIAGRAMS GENERATING THE FORM FACTOR 

We calculate the vertex loop diagrams depicted in Fig. [T], for the space-like pion momenta, 
= —Q^ < 0, and with both nucleons on-shell. 

Before renormalisation, the PTV vertex can be written as + r^(g^), which is the sum 
of the tree-level contribution from the lagrangian Eq. ([3]) and the bare (unrenormalised) 
one-loop corrections r^(g^). To render the calculated loop integrals r^(g^) ultraviolet con- 
vergent, we apply the modified minimal subtraction procedure in conjunction with dimen- 
sional regularisation [s], This procedure amounts to setting the constant A (defined 
in Appendix [X]) , which is divergent in the limit D ^ 4, to zero in the unrenormalised loops 
r^(g^). In this way we obtain ultraviolet and infrared convergent loop integrals which do 



All other conventions and definitions used throughout this paper follow Ref. 
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FIG. 1: The loop diagrams included in the calculation of the P- and T- violating ttNN form factor. 
The solid and dashed lines denote nucleons and pions, respectively. In addition, in the upper row, 
the internal dashed lines are ir or rj mesons, and in the middle row, the zigzag lines are p or u) 
mesons. The blobs represent the PTV interactions from Eq. ([3]), and the other vertices denote the 
strong interactions from Eq. ([2]). 

not depend on the regularisation parameters (see Appendix E]): the dimension D (which 
after renormalisation can be set to the physical value D = 4) and the interim cutoff A 
(which may be thought of as a scale much larger than the included hadronic masses). The 
renormalisation can be interpreted as an absorption of the infinities of the unrenormalised 
loops r'^(g^) into the bare coupling constant c^, such that the full vertex is rewritten as 
the sum of a renormalised coupling constant and a renormalised loop contribution r(g^), 
both of them free of infinities. The PTV irNN form factor F^[Q'^) is then defined by 

c^ + T{q^) = c^F„{Q^). (4) 

The sum of the renormalised loop diagrams with an intermediate pion reads 

r.(g^) = [VM') - i} + {^^(^') + 4M^C^o[M^ q\ M'; m\ M^ M']] , (5) 



where 



4M2 - g2 - g2 

m2(4M2 + g2)Co[M2, m^; m^, M^, M^] 
4M2 - g2 ' 



(6) 



with Bq and Co the two- and three-point Passarino-Veltman functions [10[, given in Ap- 
pendix |Al It can be easily verified that any dependence on the regularisation parameters A 
and A cancels out in Eq. ([5]) (as it should), thus proving the ultraviolet and infrared con- 
vergence of the renormalised loop. Using formulae from Appendix |A] and retaining only the 
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leading orders of m/M and (m/M) \n{m/M) in the Bq functions, Eq. ([5]) can be rewritten 
in a more explicit form 



+ 



+ 



647r2 167r2(4M2 - g2) 
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+ (4M^ + / dx\n\l 



'0 \ M2 ^ 

4M2(4M2 - q^)cl - m^AM^ + q^){gl + cl) Co[M\ q\ M^; m\ M^, M^] L (7) 



with Co[M2,g2 



^] given by Eq. (lASp . (The g^-independent additive constant 
is unimportant here since in presenting our results we will normalise F^{Q'^ = 0) = 1.) 
The vertex loops with the p, uj and vj mesons have a structure similar to that of the pion 
in Eq. ([5]); their contributions are given in Appendix [Bl 

Examples of the calculated PTV irNN form factor F^IQ"^) are shown in Fig. [21 where we 
display separately the contributions from the vertex loops from Fig. [T] involving vr, p, uj and 
7] mesons. The form factor depends on a particular relation among the generally unknown 
PTV meson- nucleon coupling constants Cjr, c,,, Cp and c^. This is illustrated in the four panels 
of Fig. [21 where we used PTV constants of comparable size with various sign combinations. 
The results in Fig. [2] were obtained using the pseudovector {pv) structure of the strong ttA^A^ 
vertex, as given in Eq. ([2]) . If instead we use the pseudoscalar {ps) nNN vertex to calculate 
the diagrams in Fig. [T], then only the loops containing the vector mesons p and u would be 
different from the pv case. As a result, there is an approximate relation between the PTV 
ttNN form factors obtained using the ps and the pv t^NN vertices: the form factor obtained 
using the ps irNN vertex and with a certain relation among the PTV constants c^r, c^, Cp 
and is approximately equal to that obtained using the pv vertex, but with the signs of 
Cp and Ccu changed. For example, the lower left panel in Fig. [2] shows either the form factor 



0.5c^ = — O.Sctt or (approximately) the 



0.5c„ = O.Sctt. 



calculated using the pv vertex and with Cp = c^j = 
form factor obtained using the ps vertex, but with Cp 

We can see from Fig. [2] that the contributions of the vector mesons p and u to the PTV 
ttNN form factor are larger when Cp and c^^ have opposite signs compared to c^r and of 
the pseudoscalar mesons. Furthermore, for any combination of the PTV coupling constants, 
the p and a; contributions have comparable magnitudes but opposite signs. The t] meson 

and Crj were all much 



gives the smallest contribution. If the PTV coupling constants 
smaller than c^^, then the resulting form factor would be approximated by the pion loop 
contribution alone, given by the dashed line in Fig. [2l Currently there is no strong reason to 
assume that c,r is dominant, or adopt any other particular relation among the PTV coupling 
constants (various arguments on this issue have been presented in Refs. 0, lll|). Our results 
in Fig. [2] indicate that a measurement of the PTV ttNN form factor could be very useful for 
establishing such a relation. Fig. [2] also shows that for some relations among the coupling 
constants, the form factor may have a zero at of the order ~ M"^. Although at present we 
cannot make a specific statement about the significance of such a zero, it might be important 
in some applications of the form factor. In the next two sections we will analyse in more 
detail the PTV nNN form factor in the low- and high-energy regions. 
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FIG. 2: Contributions of the meson loops to the P- and T- violating ttNN form factor F^{Q'^), 
defined in Eq. using the indicated relations among the PTV meson- nucleon coupling constants. 
The total form factor is the sum of all meson contributions and is given by the solid curves. 



IV. MEAN SQUARE RADIUS AND THE EFFECT ON THE NEUTRON EDM 

We define the mean square radius associated with the PTV ttNN interaction by analogy 
with the standard definition used for electromagnetic form factors: 
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The contributions to the radius from the loops involving particular mesons will be denoted 
by the corresponding subscripts: 

rlj) (0.19 - 0.00062 c^) fm% (9) 



6 



r|y) ^ I -0.092 ±0.015 ^ + 0.00035 c;j fm% (10) 
rjy) ^ [o.093 + 0.032— - 0.00036 c^] W, (11) 
1^) ^ ("o.OOll- 0.026 ^- 0.000018 cj] W, (12) 



where the upper (lower) sign in Eqs. (ITO!) and (fTTj) corresponds to the use of the (;js) 
structure of the strong ttNN vertex. The leading term of the pion loop contribution can be 
written analytically in a compact form 

-■W, = + = + "VM') « 0.22 W , (13) 

where we have used the Goldberger-Treiman relation Qj^Fj^ = QaM (with = 93 MeV 
being the pion decay constant and qa = 1.267 the axial coupling constant), valid in leading 
chiral order. Neglecting the (presumably suppressed) quadratic terms in Eqs. (!9l fT2l) and 
adding up the leading terms, we obtain the total mean square radius associated with the 
PTV vriVA^ form factor: 

i^rr) = i^rr). + {'rr), + {'rr)^ + {'rr), - (o-i9 ± O-Oi^^ ^ 0.032^ - 0.026^) fml 

(14) 

This shows, in particular, that if |cp,w^^|:^|c7r|, then the dominant contribution to {r'^^ 
comes almost entirely from the pion loops, in which case the explicit formula Eq. ( |T3l) can 
be used as a good approximation to the full result Eq. (HM : {r'j>j^ ~ {'rr) ' 

We can use the mean square radius to estimate the effect of the PTV irNN form factor 
on the neutron EDM D„. First we will cite here some results on Dn obtained in approaches 
similar to ours, where we will use the recent reviews 0, 12 1 and references therein. The 



upper bound from experiments with ultra-cold neutrons is D^^^ < 7.5 x 10~^^e • cm, 
where e denotes the elementary electric charge. In the Standard Model calculations, 
appears as a strongly suppressed perturbation effect: -Df^ ^ lO^^^e ■ cm. If a discrep- 
ancy between the Standard Model prediction and experiment is eventually established, it 
may indicate a P- and T-violation, or other fundamental extensions of the Standard Model. 
These possibilities are explored in many calculations wherein one uses EDM measurements 
to constrain various PTV parameters, e. g. the QCD ^-term, CP-violating quark coupling 
constants, etc. Examples of such model calculations, most relevant in the context of the 
present paper, include the following. The MIT quark bag model [l^ yields the neutron 
dipole moment D^^^ ^ 8.2 x lO'^'^^e ■ cm. The estimate from the leading chiral loga- 
rithm is D^^"'^ ^ 3.6 X IQ-^^Oe ■ cm, while the extended SU{3) chiral perturbation 
theory calculation Il5| gives to first order D^^^ (3.3 ±1.8) x 10~^^ 9 e ■ cm. In the effective 
lagrangian model j3| is evaluated in terms of the PTV meson-nucleon coupling constants 
given in Eq. P: D^^ ^ [0.14^W - 0.02 {gf^ - g^p +gf^) + Q^ 10"%°)] e ■ cm. Another 
approach based on QCD sum rules, predicts D^^ ^ 1.2 x 10~^^ 9e ■ cm. 

The aforementioned leading chiral contribution D^^°^ is obtained 0, 0, [3] from the one- 
pion loop correction to the PTV 'yNN vertex where the photon couples to the intermediate 
pion: 
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with effects of the PTV irNN form factor being neglected in Eq. (fT5|) . We assume that 
the loop integral for the neutron EDM is dominated by the region of very small four- 
momenta squared of the intermediate pions. In this region the form factor F^{Q'^) ^ 

1 — i^T^'fj^ If we set ~ m? in this expansion and substitute it into the EDM pion- 

loop correction in place of the PTV ttNN form factor, then the dependence of the form 
factor does not have to be integrated. Thus we obtain an estimate for the neutron EDM 

= Dt"' (l + , (16) 
with the correction due to the PTV txNN form factor 

^J? =» ~ {r'h) « + 0(cl, ™VM') « -2% . (17) 

where we have used Eq. ( IT3l) . We should point out that in view of the preceding simplifying 
assumptions, this is a rather crude estimate which may give only an approximate magnitude 
of the effect. A more detailed treatment must involve a careful integration of the EDM 
vertex loop including the PTV ttNN form factor (which in turn depends on the unknown 
PTV meson- nucleon coupling constants, as discussed in Section UTTl) . Such a calculation falls 
outside the scope of this paper. 



V. HIGH-ENERGY SCALE EVOLUTION OF THE PTV ttNN COUPLING CON- 
STANT FROM THE CALLAN-SYMANZIK EQUATION 

The Callan-Symanzik (or renormalisation group) equation [ol, [l^, 17 1 describes the evo- 
lution of a renormalised coupling constant Cvr(yu), defined at a typical mass scale \/—q^ ~ /i, 
where, ideally, all the masses in the problem are negligible in comparison with fi: 

= f3{c^) . (18) 

In our case it should be sufficient to take fi to be much larger than the nucleon mass M. 
The evolution is determined by the renormalisation-group (3 function /3(c7r) = jS-jr which we 
will calculate in this section. To keep the discussion transparent, we will consider the scale 
evolution of the PTV ttNN coupling constant in isolation, extracting the (3 function from 
the one-pion loops only and ignoring the evolution of the other coupling constants. Although 
these simplifications make the renormalisation group considerations somewhat formal, the 
results of this section should provide a useful approximation of the true behaviour of Cn{fi)- 
For definiteness, throughout this section we assume that C7r(/i) > 0. This choice of the 
sign does not restrict the renormalisation group considerations since, as Eq. (fT8|) shows, 
P{—Cn) = —P{cn). In general, one can consider |c7r(/i)| without any change in the following 
discussion. 

We will calculate the f3 functions for the two structures of the strong (non-PTV) tcNN 
vertex: the pv structure, given in Eq. ([2]), and the ps structure. As described in Section UTTl 
the renormalised pion loops for the PTV form factor are the same using either of the two 
structures. However, the /3 functions are different in these two will be shown below. 

We use dimensional regularisation to calculate the loop integrals at a space-time dimension D 
and apply the modified minimal subtraction procedure. In this scheme, the (3 function at one- 
loop level can be extracted as the coefficient at the term 1/(4 — D) — 7£;/2 + ln(47r)/2 = A/2 
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which enters into the bare couphng constant in order to cancel the UV divergences of the 
loop integrals [ol, 18 1. The A term appears in our calculation only through the functions 
Bq in the unrenormalised loop integrals, see Eq. (lA3p . Thus we obtain: 



(pv) 



C'jr 



in 



2 1" ' 



using the pv strong vrA^A^ vertices in the loops, and 



(ps) 



9l 



Stt' 



Stt 



2 f ' 



(19) 



(20) 



using the ps vertices. 

Integrating the Callan-Symanzik equation f|T8l) between mass scales yUi and /i2, with the 
(3 functions Eq. (fT9l) and Eq. ( !20l) . we obtain the solution 



/i2 



C.(/il)\ l + fc2(/i2) 



where 



and 



A = A^P") 



167r2 



87r2 



3.25, B 



-2.16, B 



1 



-0.013, 



-0.013, 



(21) 

(22) 
(23) 



using the pv and ps strong ttNN vertices, respectively. 

If we work at moderate mass scales /x, where \Bcl{^)/A\ <^ 1, we can retain only the 
leading order in the (3 function, 

Pn^Ac^, (24) 
for which an approximate solution of Eq. (1181) can be written 

cM ~ , (25) 

where /zq is some "initial" renormalisation point (one may choose it to be of the order of 
the nucleon mass, fio ~ M, assuming that C7r(/io ~ M) has a magnitude relevant to EDM 
measurements or to other PTV experiments). 

Consider first the calculation using the pv strong nNN vertex in the loops. In this case 
A = A(P^) > 0, see Eq. (l22l) . hence the approximate solution Eq. (125!) describes a growth of 
c^(/i) with /i. For a sufficiently large scale fi the value of c^(/x) could become so big that the 
f3 function Eq. (1191) could vanish. Then the approximate solution Eq. fl25p would no longer 
be valid. Instead, Eq. (ITSl) shows that c^(/i oo) would tend to a fixed point c* determined 
from the condition /^^^^^(c*) = 0. In the vicinity of the fixed point 



pipv) 



-Qt, ^ 16 , a = ^r-^ ^ 6.5 



(26) 



and the solution of Eq. (fTSj) describes how c^(/u) approaches c*: 



4 - cM 



(27) 
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Recalling that Cj^ is proportional to the CP-violating 6 term which could be included in a 
QCD lagrangian [2I, 0], the evolution of C7r(/i) in Eq. (!25|) is compatible with an extremely 
small ^ < 3 X 10"^" [l|, 0], as extracted from EDM measurements at low energies. In this 
scenario, obtained in our model using the pv strong ttNN vertex, an initially small value of 

grows as a positive power of the mass scale and may eventually reach a fixed value c* . 
Using Eq. with A = A(p^) ^ 3.25, /iq ^ M, c^(M) < 7 x lO^^^ ^-^ g^^^ ^^(^^^ ^ ^^/\/2 
(the point at which P^^'' of Eq. (fT9l) reaches its maximum and starts decreasing), one can 
estimate that the fixed point c* would be approached at energies above ~ 5 TeV. 

Now consider the calculation with the ps strong ttNN vertices in the loops for the PTV 
form factor. A qualitatively different behaviour of c^^lfi) obtains in this case. The jS func- 
tion Eq. fl20l) remains negative for all c^, and the general solution of the renormalisation group 
equation Eq. ( fTSi) is well approximated by Eq. (l25l) . Since A = A^^^^ < 0, see Eq. (!23|) . the 
solution Eq. ( l25l) describes a gradual decrease and eventual vanishing of c^ifJ^) as — >• 00. 

Whether one of these two contrasting scenarios of the high-energy scale evolution of 
C7r(/i) is correct should be, in principle, clarified by experiment. An important caveat here is 
that the dynamics at high energies will certainly involve hadron resonances, quark-gluon and 
other - possibly supersymmetric - degrees of freedom not included in our model (for a recent 
review, see ^ and references therein). These degrees of freedom would probably influence 
the high-energy evolution of Ct,{h), which in turn could change the discussion presented in 
this section. 

VI. SUMMARY 

We have described a calculation of the form factor associated with the parity and time- 
reversal violating (PTV) pion-nucleon coupling. The form factor is generated in our model by 
one-loop vertex corrections involving nucleon, pion, p-, uj- and //-meson degrees of freedom. 
The interaction of these hadrons is described by a phenomenological lagrangian consisting 
of a strong interaction part and a PTV part. The latter is determined by several PTV 
meson-nucleon coupling constants, all of which are in general unknown, as are their relative 
orders of magnitude. 

One result of the present calculation is that by measuring the PTV irNN form factor, 
it should in principle be possible to constrain this set of the PTV meson-nucleon coupling 
constants, since different relations among them lead to distinct dependences of the form 
factor, as shown in Fig. [2l The experimental feasibility of such a measurement is certainly 
a challenging problem in its own right, and we have left an investigation of this question 
outside the scope of the present paper. We found that the PTV nNN form factor generated 
by the loops with the strong pv nNN vertex is related in a simple way to that generated by 
the loops with the strong ps vertex: these two form factors are converted into each other by 
reversing the signs of the PTV couplings of the vector mesons. 

We have analysed the low-energy characteristics of the PTV pion-nucleon interaction by 
evaluating the mean square radius associated with its form factor. It turned out that, if 
all the PTV coupling constants are of the same order of magnitude (or at least if the pion- 
nucleon constant is not much smaller than the others), then the contributions of the p and 
u mesons almost completely cancel each other, and that of the t] meson is suppressed. As a 
result, the mean square radius can be well approximated by the pion-loop contribution alone, 
which can be given by a simple expression Eq. 0131) . This is consistent with the intuitive 
expectation that the pion cloud should dominate the low-energy features of the interaction. 



10 



The simplified estimate in Section IIVI suggests that the effect of the PTV irNN form factor 
on the neutron EDM is very small. A more definite statement can be made on the basis of 
a detailed calculation of the loop integrals for the EDM, including the full Q^- dependence 
of the PTV ttNN form factor, which was not pursued in the present paper. 

It has been shown elsewhere [g], that the PTV pion-nucleon coupling constant is 
proportional to the CP-violating 6 term which can be added to the QCD lagrangian. The 
current estimates of 6 based on measurements of EDMs of the neutron and various atoms 
and molecules, yield an extremely small upper limit \6\ < 3 x 10"^", which presents a 
long-standing puzzle [2|. Since the value of 6 may be energy-scale dependent (for example, 
yielding appreciable CP- violating effects at sufficiently high energies), we have considered 
in our model possible scenarios of the evolution of the PTV nNN coupling constant c-jrifi) 
with the energy scale fi. To this end, we employed renormalisation group techniques, with 
several simplifying assumptions. Being used in the framework of a model, results of such an 
analysis are to some extent speculative. Nevertheless, we found in Section |V] an intriguing 
possibility that, in the calculation with the pv strong irNN vertex, the PTV irNN coupling 
constant relevant to very high energy scales may be significantly larger than the tiny value 
extracted from the low-energy EDM measurements. In the ultraviolet limit /x — > cxd, Ct^Ih) 
may eventually reach a fixed value c*, Eq. fl26l) . of a "natural" size. 



APPENDIX A: PASSARINO-VELTMAN FUNCTIONS 

In this appendix we give the scalar Passarino-Veltman functions 10, To'l Bq and Co enter- 
ing into the expressions for the loop integrals. The functions are defined as D-dimensional 
integrals 



4-Z) 



Bo\p';ml,ml]= ^^""^^^ ^ f d"" q {[q^ - + tO][{q + pf - ml + tO]} ' , (Al) 

Co[pl,pl,{pi+P2f;ml,ml,ml] = 

^^""^^a" / ^""l { il^ - + [(^ + Pi)' - + «0] [{q + pi+ P2f -nil + zO] , (A2) 

where A is an arbitrary cutoff mass introduced to keep the dimensions of and Co inde- 
pendent of D (upon renormalisation, the loop integrals do not depend on A, and the limit 
D — i> 4 can also be taken without encountering divergences). Eqs. (lAll) and (]A2I) can be 
reduced to one- and two-dimensional integrals which are either evaluated numerically or 
expressed via logarithm and dilogarithm functions. We have 

Bo[p^ m\, m^] = A + Bl\p'- m\, m^] , (A3) 

where 

2 

A = ^ — — - 7ij + ln(47r), 7^; ^ 0.5772 (Euler's constant). (A4) 

is a divergent (in the limit D ^ 4) constant typical for the modified minimal subtraction 
procedure, and 

Bl\v'-m\,m'^ = -J^ dx\n i-^ 1^ \ + 0{D - A) , (A5) 
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is a finite (ultraviolet convergent) one-dimensional integral. 

In our loop integrals we encounter only Co functions with pf = {pi + ^2)^ = M"^ and 
p\ = q^i which can be reduced to convergent two-dimensional integrals 

Co[M\q\M^-mlmlm^^ = 



—J dx J dy |m^(x + uY — xyq^ — M'^{x + y) — x{m\ — m^) — y{m\ — mg) + — ioj . 

(A6) 

The Bq and Co functions from the pion loops Eq. ([5]) can be further simplified: 



, M /•! / q^x(x - 1) ^ 
-21n— - / dx\n[l + ^—— — '--iO 
A Jo I M2 



(A7) 



M2(l-x2) 



\n{xg) 



Co[M\q\M';m\M\M' 

1 Tfl 

21n(l + x,)--lna;,-21n— 



+ — + 2Li2(-x, 



9/ f ' 



(A8) 



where x„ 



(g2 _ 4M2)/(g2 + iO) - 1 ■ J(g2 - 4M2)/(g2 + iO) + 1 and the diloga 



-1 



rithm (Spence) function is 
Li2(x) 



1 ln(l — xy) 

dy 

y 



larefl — x)\ < tt ) . 



The following special cases are useful in deriving Eqs. (ITj) and (|T3 



Bl[q^-M\M'' 



, M g2 
-21n— + ^ 



A 6M2 



M4 



Bl\M^\m^,M^] = -2 n— + 2 1 arctg— ^ n— + 0\ — 



m 



m 



m 



2 \ 



Co[M^ 0, M^; m\ M\ = ^,^^^^0\ — \. 



(A9) 



(AlO) 

(All) 
(A12) 



APPENDIX B: LOOP INTEGRALS WITH THE p, uj AND r] MESONS 

In this appendix we give explicit formulae for the renormalised loop diagrams with p, to 
and ?7 mesons, where for brevity we will retain only the dominant linear terms ~ Cj^^p^uj^-q-, 
marking them with a superscript "(1)"- The contribution from the 77 meson equals 



r«(g^ 

where 



2 2 

C-qg-Kg-qQ („2\ D r/\/f2. /\#2 _2ll , ^T^dr] 

'l67r2 



(Bl) 



^7r2(4M2 - g2 

V^{q') = B,[q'- M\ M^] + mj Co[M\ q\ M'; m% M^ M^] . 



(B2) 
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The contributions from the u and p mesons read 

+ 2c^g^q^ (8M2 + 12k^M^ - 2K^ml - 3K^q^) } Bo[q^; M\ M^] 

+ 64M%^(4M^-g^) {''^^'^ (^^^' + ^^^'^^^ + 96«.M^m5 

+ 324 M%i2 + 4/^2 ^fVif, - 48M^g^ - 24Ka;M2m^g2 - IG^M^m^g^ 

+ '^>^?' + SM^g^ + 2k>^^^) 

+ 2c^g^mlq^ (aM^ + 8«:,M=^ - «:,m^ - 2K^q^) }Co[M\ q\ M^; m^, M^, M^], (B3) 



+ 32M%^ |l + ^o[0,m ,M ]| ^^^^ Bo[0,m 

. g/mCpg-K \ TIT? . 2 2i d Fn /\/r2 n#2il , QpnCpg^iM —m) „ 2l 

" 32M%2f4M2 _ ^2) {c-^P?' ('^P^' - - 4«P^' - «p^?) 

+ 2cpg^ {2m^M'^ + 2M2m2 - m^q^ - 2M'^q^) } Bq[M^- m^ M^] 

+ 32M%2flM2 - q-) {'^^""^^ (^^' + ^''^^^ + '^^^^ - '^''^') 

+ 2cpg^ {2m^M^ + 2M2mJ + 2M^(^ - my) } Bo[M^; M^, 

~ 47r2(4M2 _ g2) {c-^P (4^' + '^P^p + 5' + «^P5' - '^P^') 

+ Cpg^ {m'^ - "^p + } -Bo[?^; 

+ g^2(4Af2 _ g2) (SM^mJ - 8m^M^ + K^m^ - K^m^ + Avr?q^ 

+ 2Kpm}q^ - 8M^g2 _ ^^^4j 

+ Cpg, (m^ + - g^) (m^ - + g^) } Co[M^ g^ M^; M^ m^ , (B4) 
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with Bq and Cq defined in Appendix [X] To facilitate tlie calculation of the loop integrals, 
we used the computer package "FeynCalc" (20| . 
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